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Abstract 

This paper considers the control landscape of quantum transitions 
in multi-qubit systems driven by unitary transformations with single- 
qubit interaction terms. The two-qubit case is fully analyzed to reveal 
the features of the landscape including the nature of the absolute 
maximum and minimum, the saddle points, and the absence of traps. 
The results permit calculating the Schmidt state starting from an arbi- 
trary two-qubit state following the local gradient flow. The analysis of 
multi-qubit systems is more challenging, but the generalized Schmidt 
states also may be located by following the local gradient flow. Finally, 
we show the relation between the generalized Schmidt states and the 
entanglement measure based on the Bures distance. 
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1 Introduction 



The topology of quantum control landscapes is important because it estab- 
lishes the general features of the control behavior generated by applying 
external fields [HE]- The landscape for quantum transitions, assuming com- 
plete controllability, was analysed with the conclusion that there are no traps 
HI [5j [6] that could hinder achieving the highest possible control outcome. 
This paper studies the problem of describing the landscape of quantum tran- 
sitions driven by local unitary operators, i.e., those acting on one qubit at a 
time, for multi-qubit systems [7]. 

The Schmidt states, defined for pure bi-partite systems, are important 
because of the insight they can provide about entanglement. The Schmidt 
states were generalized in [HI [9], in order to treat multipartite systems. This 
paper will show how to obtain the canonical form of the generalized Schmidt 
states by following the local gradient flow. This technique ultimately leads to 
a method to measure the entanglement of pure systems based on the optimal 
implementation of local unitary operations as a subset of the more general 
classical operations and classical communication protocols as it was pursued 
with other methods [TQl [TTJ. 

It is convenient to define the following bracket operation 

Wo = ^Tr[X + Xt] (1) 

The comparative fidelity between two density matrices, when at least one of 
them is pure is (poPt)o- If one of the states is driven by a unitary operator, 
then the cost function can be written as 

F = (tfp Up T ) , (2) 

This expression has the same form as the cost function for the optimization 
of the expectation value of an observable O [5], 

Jt = (Up tfO) , (3) 

which was the subject of prior landscape studies [2J. The fidelity function 
for the state transfer can be rewritten as 

F = (poUp T U^) . (4) 
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An infinitesimal transformation of the unitary operator can be expressed as 

£/_> U' = Ue 5A = U(1 + 6A), (5) 

with 5 A being an anti-Hermitian element that lies in the corresponding Lie 
algebra, so that an infinitesimal variation of U becomes 

5U = U5A, (6) 

which can be used to calculate the first order variation of the fidelity as 

6F=(p U[6A,pT]rt) . (7) 

A subsequent manipulation results in 

5F = (lp T , tfp U]6A) = (lp T , tfp U)tf6U) (8) 

thereby identifying the gradient as 

Gradi = U[U* pall, pr], (9) 

with the corresponding gradient flow equation 

^L = u[u^ Po u,p T ]. (io) 

as 

The fidelity can be expanded up to second order to obtain the quadratic 
form for the Hessian 

5 2 F = ({pr, U^ Po U}(5A) 2 ) - 2(lP Po U6Ap T 6A) , (11) 

where {, } stands for the anti-commutator. This quadratic form is simplified 
at the critical points where the gradient (I10p is zero 

6 2 F\ C = 2((p T U^ Po U(SAf) - (U^p USA PT 6A) Q ). (12) 

The local gradient flow is found by eliminating multi-qubit terms in SA, 
such as <7 3 ®0"3 and leaving single qubit terms, such as <r 3 ® l 2x2 or 1-2x2 ®°3- 
In this way, only strictly localized interactions are involved as happens in 
classical mechanics. Defining V as the projector that eliminates multi-qubit 
terms, the variation of the unitary operator with the corresponding local flow 
is 

5U = UV5A. (13) 
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The projector V is easily calculated by tracing one-qubit terms. For example, 
the two-qubit projector is 

1 3 

V — - ^2 T r [ ■ °o <8> &j]vo <8> Oj + Tr[ ■ a i <g> a ]aj <g> a , (14) 

3=1 

with (T = 1-2x2, so that V8A is constrained to the six-dimensional Lie algebra 
su(2) x su(2) C su(4). The first order variation subject to the local flow 
becomes 

5F = ([p T , U^p U]VU^5U) = (V (Ipt, U^ Po U}) U^5U} , (15) 
which results in the following local gradient 

Grad[ ocal = UV[U^p U, p T \. (16) 



2 Two-Qubit Systems 

The Schmidt states play an important role in the quantification of the entan- 
glement of two-qubit systems. We will show their importance in describing 
the quantum landscape characterized by the local gradient flow and then 
calculate the Schmidt state of a given entangled state by following the local 
gradient flow (excepting the maximally entangled state). 

Consider the landscape where the target state is a Schmidt state denoted 
as pt = Ps(Q)- The Schmidt states for two-qubit systems can be parametrized 
with a single variable as 

|^ s ) = cos(^/2)|TT) + sin(^/2)|||), (17) 

whose corresponding density matrix reads 

/cos 2 (#/2) isintf \ 



Ps(0) 





y \ sin^ sin 2 (#/2)y/ 



(18) 



with < 9 < 7T, in the standard basis {| |T ), I U ), I IT ), I II >}• 
The critical states p c = U\p§U c obey the following equation 

V[ Pc ,ps(9)} = (19) 

It can be shown that this equation is satisfied by critical states that fall into 
one of the following two cases 
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Another Schmidt state p c = ps{<t>)- 111 this case, the eigenvalues of the 
Hessian around the critical points are either negative or mixed, with 
the following explicit form 



MM) 



o 



\ 



— 1 — cos(# — (f)) — sin 9 — sin <; 

— 1 — cos(# — 0) — sin 9 — sin <; 

— 4 sin # sin 
—1 — cos(# — 4>) + sin 9 + sin <; 
y — 1 — cos(# — <fi) + sin 9 + sin <f> J 



(20) 



For each critical state with a negative spectrum h(6o,<j)o), there is an- 
other one with a mixed spectrum h(9o,7c — <f>o). Conversely, for each 
critical state with a mixed spectrum h(9 , <f> ), there is another one with 
a negative spectrum h(9o,Tt — O )- So, for each initial state there is a 
pair of critical states that can be reached by following the local gradient 
flow, such that one of them is a saddle point and the other is an stable 
maximal point. If the initial state is separable, the two possible critical 
states are given by ps(0) or ps{^)- 

The critical sub-manifold spanned by the basis {| T! )> I IT )} with the 
following explicit form of the critical state 



(21) 



(22) 



p c = x\ |t)<IT| + (i-z)|TI}(TI 

where the eigenvalues of the Hessian are 

(l - -2x) 2 cos 2 # + sin 2 e\ 
1 - y/(l - 2x) 2 cos 2 9 + sin 2 9 


1 + v/(l-2x) 2 cos 2 ^ + sin 2 ^ 
1 + v/(l-2a;) 2 cos 2 ^ + sin 2 ^ 
V o / 

which corresponds to a positive spectrum, associated with the mini- 
mum. 

Based on the features of the critical points we can state the following 
theorem 
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Theorem 1 The fidelity landscape between a pure separable state p and a 
target Schmidt state ps(0) (with 9 ^ it/ 2) has saddle points but no traps. 
Moreover, the separable states that maximize the fidelity converge to either 
I TT ) or I II ) depending on the target state as they follow the local gradient 
flow, according to the following formula 

, lmC /V„t/=( ITT><TTI °<"<^ 2 , (23) 

This theorem is a direct result of the fact that those limiting states are the 
only Schmidt states with zero entanglement. Moreover, we can also say that 



Corollary 1 For pure states, the maximum fidelity between an entangled 
state and a separable state can be calculated from the corresponding Schmidt 
state \ip s ) = cos(0/2)| TT ) + sin(0/2)| || ) as 



cos 2 (0/2) 9 < tt/2 
sin 2 (0/2) tt/2 < 9 < vr 



T{9) = max F = { J t 7 n . _ ( 24 ) 



The maximum fidelity F{9) can be used to calculate the Bures distance 
as the entanglement measure, which satisfies all the features required for a 
good entanglement monotone [12], [13] . In the present case of pure two-qubit 
systems the entanglement formula is 

E B {p) = 2 (l - v 7 ^) • (25) 

As a first example, Figured] shows the fidelity of the states following the local 
gradient flow for the initial separable state described by 

p = e £"o®"i| TT)(TT |e-^ CT0 ® CTl (26) 

with ps(7r/4) as the target state and | TT ) as the limiting state. The next 
example considers the following entangled initial state 

p = e ^ CT2 ® CT0 e ^ <72 ®°' 2 | TT)(TT |e~i^ <T20CT2 e~^74 CT2 ® <TO (27) 

driven by the local unitary flow with ps(tc/4) as the target state, and the 
following limiting Schmidt state 

/0. 793893 0.404508\ 





\0.404508 0.206107/ 



lim tfpoU 



(28) 
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Figure 1: Fidelity curve of the states following the local gradient flow for 
the initial separable state (121)1) with ps(7r/4) as the target state. The fidelity 
never reaches 1 but attains the global maximum associated with the limiting 
state | tt )• 

Almost any Schmidt state can be used as the target state in order to drive 
the local gradient flow, excepting those with 9 = {0, 7r/2,7r}, because of 
convergence issues. For example, Figure [2] shows how the arbitrary state 
(127)1 approaches its Schmidt state for the range of target Schmidt states. 

The local gradient flow was driven by employing target Schmidt states, 
but the landscape is invariant under the application of local unitary opera- 
tions on both the initial and target state. The local unitary transformations 
include local phases, that are able to change the phase of the Schmidt states. 
This means that the general stable critical states are Schmidt states with 
the possibility of extra phases. For example, consider the following arbitrary 
entangled state made from a Schmidt state and local unitary transformations 

p E = e f ^^ogf-o^i^^^^-f <xo^i e -f- 2 ^o, (29) 

The initial separable state is taken as pi = | i t )(TT I- The local gradient flow 
converges to a separable unitary operator U c with the following corresponding 
separable state 

(30) 



- s 
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Figure 2: Fidelity of a random entangled state moving towards its Schmidt 
state as a function of the target Schmidt state ps(0) employed to drive the 
local gradient flow. The dashed line represents the fidelity of the initial 
random state with respect to its Schmidt state and each subsequent curve 
corresponds to another step in the approach by following the local gradient 
flow. The figure suggests that the arbitrary state never reaches its corre- 
sponding Schmidt state when the gradient employs the target states ps , (0), 
P 5 (tt/2) and p<?(tt). 




8 



This state can be diagonalized by the following local unitary operator 



T 



l/y/2 1/V2\ 
-l/y/2 l/y/2) 



i/y/2 i/y/2\ 
-l/y/2 l/y/2) 



(31) 



such that T^p c T = | tt )(TT I- This suggests that T could be used to reduce 
Pe to its expected Schmidt state ps(7r/4), but instead we obtain a Schmidt 
state with an extra phase —% 



( 



T^p E T 



cos 2 (f) 








V cos (I) sin (!) 










-zcos(f)sin(f)\ 





sm 



(32) 



(I) I 



This extra phase can be eliminated by the use of local phase transformations, 
which otherwise leave the absolute value of the components of the density 
matrix invariant. 



3 Three or More Qubit Systems 

The entanglement in a two-qubit system can be minimally characterized by 
a single variable as shown in the Schmidt state. The number of variables 
needed to parametrize a n-qubit system is 2 n+1 — 2 up to a global phase, 
and the number of variables to parametrize a single qubit is 3n, thus, the 
minimum number of variables needed to parametrize the entanglement of an 
n-qubit system is 

N E = 2 n+1 - 2 - 3n, (33) 

which is five for three-qubit systems. The canonical form of the generalized 
Schmidt states is important because of the information that can be obtained 
about entanglement [Til [T5| [16j [TJ]. A canonical form of the generalized 
Schmidt state for three-qubits was introduced in [8] as 

\1>s) = Ail TTT) + A 2 e«*| Til) + A 3 | IT!) + A 4 | lit) + A 5 | III), (34) 

with Aj > 0, <p > and |Aj| 2 = 1. The canonical form of the generalized 
Schmidt state for n-qubit systems was given in [9j indicating that the missing 
basis elements in the generalized Schmidt state are 

I ITT ••• T), I TIT ••• T), I TTI ••• T),-l TTT ••• I)- (35) 
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The landscape of multi-qubit systems is richer and more complex than the 
two-qubit case. Considering the case where the initial state is separable and 
following the reasoning in [9], we can always demand that Ai > Afc. However, 
the analysis is simpler if we relax some generality and demand that Ai > Afc, 
for k > 1. The variation of the fidelity can be written as 



6F = 6(y\i> s ) = 2Re[(V\ip s )(ip s \6\V)] 



(36) 



The canonical form in (j34j) indicates that if we start with a generic separable 
state |^) = \ipi) ® \ip2) <8> 1^3) an d allow local unitary transformations, the 
isolated maximum fidelity is achieved at the critical state |\& c ) = | Ttt }■ The 
first order variation under local unitary transformations is made of a linear 
combination of basis elements with at most one qubit reversed, 



8\ TTT) = (l + i6i)\ TTT) + S 2 \ ITT) + fc| TIT) + S 5 \ TU), 



(37) 



with 81 G R, 62 G C, 83 G C, 85 G C. We can use this variation in order 
to evaluate 8F given by (|36|) at the critical state |^ c ) and verify that it 
is a stationary point, thus justifying the canonical form of the generalized 
Schmidt state. The missing basis elements (1351) form a critical sub-manifold 
associated with the fidelity minimum of zero value. The generic identification 
of the remaining critical states is difficult and depends on the specific Xj 
values. However, if Xj > 0, then there are no additional critical states because 
the aforementioned critical states exhaust all the possibilities to obtain 5F = 
0. 

As a concrete example, consider calculating the generalized Schmidt state 
of the following arbitrary state 

/0.3 + 0.H\ 
0.2 
0.3 
0.3 
0.4 
0.2 
0.5 



|V>T> 



(38) 



0.77/ 

Following the same procedure used in the two-qubit case, we use the local 
gradient flow to calculate the optimized separable state \i/) c ) that maximizes 
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the fidelity K^tIV'c)! 2 , starting form an initial separable state (e.g. | T T T ) ) - 
The optimized state \ip c ) can be diagonalized using a local unitary transfor- 
mation. Applying the same local unitary transformation to the target state 
|'0t)('0t| we obtain 



/ 



1^5) 



0.986657 



-0.125609 - 0.0245643i 


0.0151643 - 0.0312796i 
0.0703562 + 0.0477398i 
\-0.0138602 + 0.038707H/ 



(39) 



(40) 



which is almost in the canonical form (134"|) . The first component can always 
be put in real form by choosing a suitable global phase. The remaining 
procedure is to employ the three available local phase transformations in 
order to eliminate the phase of last three components to finally obtain 

/ 0.986657 \ 




-0.125609 - 0.0245643i 


0.0347616 
0.085024 
\ 0.0411138 J 

which we ascertain to be the global maximum because |Ai| is greater than 
the rest of the components. The local phase transformations do not change 
the absolute value of the components of the column spinor, so, it is easy to 
verify that, for example, in the last component | — 0.0138602 + 0.0387071z| = 
0.0411138. 

The procedure to calculate the Schmidt state can be used to calculate the 
Bures distance as an entanglement measure if |Ai| is greater than the rest of 
the components. In this case the formula of the Bures distance as a measure 
of entanglement is simply 



£ B (p) = 2(l-|Ai| 



(41) 
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The study of higher multi-qubit states follows along the same general lines of 
the three-qubit state. Thus, we are able to calculate the generalized Schmidt 
state as well as the Bures distance as a measure of entanglement for most of 
the cases where Ai results in a value greater than the rest of the components. 

4 Conclusions 

The landscape of local quantum transitions for two-qubit systems is well 
suited for optimization through the gradient flow because of the lack of 
traps. We showed how to extend these results to muli-qubit systems and 
presented an example on how to calculate the generalized Schmidt state for 
three-qubits. The local gradient flow can be easily applied to higher multi- 
qubit systems and even though we could not give a complete analysis of 
the landscape, a criteria was presented to establish if the global maximum 
was attained. A generalization of this analysis to mixed multi-qubits is de- 
sireable, but this is a much more challenging problem because of the severe 
limitations that unitary transformations present. 
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